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Abstract 

A unique probabilistic theory is described to predict the properties of nanocomposites. The simulation 
is based on composite micromechanics with progressive substructuring down to a nanoscale slice of a 
nanofiber where all the governing equations are formulated. These equations have been programmed in a 
computer code. That computer code is used to simulate uniaxial strengths properties of a mononanofiber 
laminate. The results are presented graphically and discussed with respect to their practical significance. 
These results show smooth distributions. 


Introduction 

The research in the nanoscale technology has exploded over the recent past. An indication of this 
explosion is that the SAMPE (Society of Aerospace Material and Processing Engineers) Conference is 
devoting four sessions of about six papers each in the last 3 years. These papers cover practically all 
current research activities. The majority of the research is devoted to processing because of the difficulties 
involved in making a useful material (ref. 1). A few investigators have been fortunate to make some 
testing samples, which they subsequently tested to obtain limited data (ref. 2). A few other investigators 
researched the characterization of fatigue (ref. 3) and creep (ref. 4). A couple of papers explored the 
construction of nanocomposites for rocket ablative material (ref. 5) and for carbon nanotubes for adaptive 
structures (ref. 6). One paper ventured to describe a computer simulation of macroscopic properties of 
carbon nanotubes polymer composites (ref. 7). However, there are no results of what special macroscopic 
properties are included. Reference 7 shows one stress strain curve and citation of several references. One 
recent article (ref. 8) describes multiscale modeling and simulation of nanostructural materials from 
atomistic to micromechanics. This article does not include information on nanocomposites, but it 
mentions that mechanistic models will be needed in the end. It is becoming abundantly clear that no 
holistic approach has been used to investigate the mechanistic prediction of all nanocomposite uniaxial 
properties. 

In this paper a unique mechanistic method is described to probabilistically simulate five uniaxial 
strengths and the transverse modulus of a mono nano fiber uniaxial composite. The mechanistic 
deterministic simulation of all uniaxial properties is described in reference 9. 

Fundamentals 

The fiber alignment with uniform dispersion is not met in nanocomposites. It is assumed herein that 
the fibers are aligned only for predicting “point” through-the-thickness properties. The fussiness can be 
simulated by estimating the angle of single fibers through the thickness. Therefore, it is assumed that an 
aligned unidirectional typical section of a nanocomposite is as illustrated schematically in figure 1 on the 
left 1(a). A nanoply is schematically shown in figure 1 on the right 1(b). It is interesting to note that in 
substructuring into slices, the monofiber nanoply is not constrained by the maximum fiber volume ratio, 
even though the monofiber was assumed to be in a square array with a limiting fiber volume ratio of 
about 0.78. The input includes the constituent material properties, tables I and II, the fabrication 
parameters, environmental, and the loading conditions. 
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(a) (b) 

Figure 1. — Unidirectional nanocomposite typical section, (a) Nanocomposite, (b) Nanoply. 


TABLE I. — T300 GRAPHITE NANOFIBER (PYROGRAF II) PROPERTIES 
[Conversion Factors: 1 10 nm = 2.756xl0 -6 in.; psi = 6.89 Pa; lb/in. 3 = 1 146 kg/cm 3 ; 
in./in./°f = (2/5); cm/cm/°c; btu = 1055 joules.] 


Description 

Symbol 

Value 

Units 

Number of fibers per end 

Nf 

1.0 

number 

Filament equivalent diameter 

df 

2.756xl(T 6 

in. 

Weight density 

Rhof 

0.064 

lb/in.* *3 

Normal moduli (11) 

Efl 1 

l.OxlO 9 

psi 

Normal moduli (22) 

Ef22 

7.0xl0 7 

psi 

Poisson’s ratio (12) 

Nufl2 

0.2 

Nondimensional 

Poisson’s ratio (23) 

Nuf23 

0.25 

Nondimensional 

Shear moduli (12) 

Gfl2 

5.0xl0 7 

psi 

Shear moduli (23) 

Gf23 

3.5xl0 7 

psi 

Thermal expansion coefficient (11) 

Alfafll 

-5.5xl0“ 7 

in./in./°F 

Thermal expansion coefficient (22) 

Alfaf22 

5.6xl0“ 6 

in./in./°F 

Heat conductivity (11) 

Kfl 1 

444.0 

Btu/hr/ft 2 /°F/in. 

Heat conductivity (22) 

Kf22 

4.0 

Btu/hr/ft 2 /°F/in. 

Heat capacity 

Cf 

0.22 

Btu/lb/°F 

Dielectric strength (11) 

KeFll 

0.0 

V/in. 

Dielectric strength (22) 

Kef22 

0.0 

V/in. 

Dielectric constant (11) 

Gamma 1 1 

0.0 

in./V 

Dielectric constant (22) 

Gamma22 

0.0 

in./V 

Capacitance 

Cef 

0.0 

V 

Resistivity 

Ref 

0.0 

Q-in. 

Tensile strength 

SfT 

8.0x10 s 

psi 

Compressive strength 

SiC 

6.0x10 s 

psi 

Shear strength 

SfS 

4.0x10 s 

psi 

Normal damping capacity (11) 

psil If 

0.38 

%Energy 

Normal damping capacity (22) 

psi22f 

6.3 

%Energy 

Shear damping capacity (12) 

psil2f 

3.34 

%Energy 

Shear damping capacity (23) 

psi23f 

6.3 

%Energy 

Melting temperature 

TMf 

6000.0 

°F 
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TABLE II.— INTERMEDIATE MODULUS HIGH-STRENGTH MATRIX (EPOXY) 
[Conversion Factors: 110 nm = 2.756x1 0“ 6 in.; psi = 6.89 Pa; lb/in 3 = 1 146 kg/cm 3 ; 


in./in./°F = (2/5); cm/cm/°C; btu = 1055 joules] 


Description 

Symbol 

Value 

Units 

Weight density 

Rhom 

0.044 

lb/in.* *3 

Normal modulus 

Em 

500000.0 

psi 

Poisson’s ratio 

Num 

0.35 

Nondimensional 

Thermal expansion coefficient 

Alfa m 

3.6xl(T 5 

in./in./°F 

Heat conductivity 

Km 

0.008681 

Btu/hr/ft 2 /°F/in. 

Heat capacity 

Cm 

0.25 

Btu/lb/°F 

Dielectric strength 

Kem 

0.0 

V/in. 

Dielectric constant 

Gammam 

0.0 

in./V 

Capacitance 

Cem 

0.0 

V 

Resistivity 

Rem 

0.0 

Q-in. 

Moisture expansion coefficient 

Betam 

0.0033 

in./in./%moisture 

Diffusivity 

Dm 

2.16xl(T 7 

in.** 2/hr 

Saturation 

Mm 

0.0 

%moisture 

Tensile strength 

SmT 

15000.0 

psi 

Compressive strength 

SmC 

35000.0 

psi 

Shear strength 

SmS 

13000.0 

psi 

Allowable tensile strain 

eps mT 

0.02 

in./in. 

Allowable compression strain 

eps mC 

0.05 

in./in. 

Allowable shear strain 

eps mS 

0.035 

in./in. 

Allowable torsional strain 

eps mTOR 

0.035 

in./in. 

Normal damping capacity 

psiNM 

6.6 

%energy 

Shear damping capacity 

psiSm 

6.9 

%energy 

Void heat conductivity 

Kv 

0.0012 

Btu/hr/in./°F 

Glass transition temperature 

Tgdr 

420.0 

°F 

Melting temperature 

TMm 

0.0 

°F 


The properties prediction is expedited by the following geometric diagrams: An exploded view of 
nanoscale isolation of a typical part is shown in figure 2 with nanoscale dimensions. A single nanofiber 
schematic with substructuring is shown in figure 3(a), and a typical subslice is shown in figure 3(b). 

A nanosubply with its corresponding stresses is shown in figure 4. The nanomechanics predictive 
equations are derived by using figure 4. The equations used are all programmed in ICAN/JAVA (ref. 10). 
A limited set of equations used for this paper are summarized in the appendix. 

Prior to describing the results obtained, it is instructive to describe the interphase and how it is 
modeled. The schematics in figure 5 show a vertical section, upper figure part, with unit thickness of the 
nanocomposite and a single fiber in it. As can be seen in the slice, lower figure, the fiber interphase is 
represented by a series of progressively larger volume voids starting with the smallest near the matrix 
interface and ending with the largest in the fiber interface. It can be visualized that the stress in the matrix 
will be magnified because of the voids. This magnification is shown in figure 6 for a specific 
nanocomposite with 0.05 fiber volume ratio and with void volume ratio varying from 0.05 to 0.4. The 
interesting point to note in the lower part of figure 5 is that the matrix is continuous even though it is 
filled with progressively larger voids; otherwise the stresses will not be continuous in the matrix. 

It is instructive to elaborate a bit further with the geometry of figure 5, lower part. In order to fill up a 
conventional ply of 0.005 in. thick and a width of 1 in., it will require about lxlO 6 fibers, a very large 
number indeed. The magnification factor of the voids effect in the interphase is show in figure 6. As can 
be seen in figure 6, the magnification factor increases from a value of about 1.1 to a maximum of about 2. 
Therefore, the maximum void effect will be nearest to the fiber interface. 
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■ Nanoscale 

diameter 
(2.756x10^6) 

Figure 2. — Nanoscale isolation of a typical part (units are in in.). 



Figure 3. — Nanofiber substructuring, (a) Several slices through the thickness, (b) Nanofiber sliced. 



Figure 4. — Nanostresses on a nanosubply (units are in in.). 
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Nanofiber in nanomatrix; dimensions in nanoinches 


Nanofiber 
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Figure 5. — Vertical section a composite nanocell through nanofiber center. 
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Results and Discussion 


In this section the probabilistic results are presented and discussed starting with the large voids in the 
interphase. The probabilistic void magnification factor is shown graphically in figure 7. It can be seen in 
figure 7 that the larger the void content the greater the deviation. The left most figure is closest to the 
matrix interphase interface while the right most curve is closest to the interphase interface. The respective 
scatter is about 0.1 for the curve closest to the matrix to about 1 for the curve closest to the fiber. The 
corresponding sensitivities are shown in figure 8. It can be seen in this figure that the void sensitivities on 
the magnification factor is large. The probabilistic void effects on the uniaxial strengths are plotted in 
figure 9. Figure 9(a) shows the spread in the longitudinal tensile strength; figure 9(b), in the longitudinal 
compressive strength; figure 9(c), in the transverse tensile strength and figure 9(d) in the transverse 
compressive strength. It can be seen in figure 9 that the distribution for the two longitudinal strengths is 
relatively large. It is from 150 to 650 ksi, for tensile strength and with a distribution of about 500 ksi, and 
for the compressive strength is from 140 to 500 ksi or a distribution of about 360 ksi. The corresponding 
probabilistic sensitivities are plotted in figure 10 for tensile and figure 1 1 for compressive. It can be seen 
in these two figures that there is no difference in the sensitivities for the three probabilities. 

The probabilistic intralaminar shear strength is plotted in figure 12. The distribution in this strength 
is from about 6,000 to -16,000 psi or -10 ksi spread. It is a relatively wide distribution from lowest 
probability to the highest. The corresponding probability sensitivities are plotted in figure 13 for uniaxial 
nano transverse tensile strength. Note that these probabilistics are for 0.0001, 0.50, and 0.9999. They are 
about the same and may be easily interchangeable as well as for three fiber volume ratios. 

The corresponding sensitivities for the nano uniaxial compressive strength are plotted in figure 14 for 
three probabilities, as was the case for the transverse tensile; these sensitivities are also the same and can 
easily be interchangeable. The corresponding sensitivities for the nano intralaminar sheer strength are 
plotted in figure 15. As was the case for the previous sensitivities, these sensitivities are also about the 
same and can easily be interchanged. The probabilistic transverse modulus is plotted in figure 16. The 
scatter is from 225 to -650 ksi or a distribution of -425 ksi. 
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Figure 8. — Voids sensitivities on the interphase magnification factor. 



Figure 9. — Probabilistically plotted nano uniaxial strengths. 
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Figure 10. — -Probabilistic sensitivities for nano longitudinal uniaxial strength. 
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Figure 11. — Probabilistic sensitivities for the nano compressive uniaxial strength. 
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Figure 12. — -Probabilistically plotted intralaminar uniaxial shear strength. 
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Figure 13. — Probabilistic sensitivities for nano uniaxial transverse tensile strength for three different probabilities. 


NASA/TM— 2007-2 14847 


9 



Probability Levels 



□ sm:'; 



OI-K7? 

ague 



Figure 14. — Probabilistic sensitivities for nano uniaxial transverse compressive strength for three different probabilities. 
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Figure 15. — Probabilistic sensitivities for nano intralaminar shear strength for three different probabilities. 


■«r 

■ hM 
P«LIM’ 
□ Gris 

■ 5M£ 


NASA/TM— 2007-214847 


11 



Figure 16. — Probabilistically plotted nano transverse modulus. 


Concluding Remarks 

The salient remarks from an investigation to characterize an aligned mono fiber nanolaminate are as 
follows: 

1 . The characterization for the nanolaminate (composite) was based on a series of progressive 
substructuring down a sliced single-diameter fiber. 

2. The theoretical development and all the equations are included in a computer code called 
ICAN/JAVA. 

3. The characterization includes two fabrication parameters, 5-nano uniaxial strengths and the 
transverse modulus. 

4. The nanolaminate investigated consists of single nanofiber laminate with 0.05 fiber volume ratio. 

5. The effects of the interphase are especially important and are represented by progressively large 
amounts of voids from the matrix interface to the fiber interphase. 

6. The probabilistic evaluation characterizes the effects of uncertainties in all participating variables. 

7. The voids uncertainties indicate as the void volume ratio increases the distribution increases as 
well. 

8. The voids contribute significantly to matrix dominated strengths. 
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Appendix — Equations Used in the Nanomechanic Characterization 


Partial volumes 

Ply density 
Resin volume ratio 

Fiber volume ratio 
Weight ratios 
Ply thickness (S.A.) 

Interply thickness 

Interfiber spacing (S.A.) 
Contiguous fibers (S.A.) 


k f + k m + k v — 1 


^Void (k v ) 



Matrix (m) - 1 


Fiber (k f ) 


P t — k f P / + k m p m 

k m ={ 1 - k v )/[l + (p m /p/)(lA m - 1)] 

kf={l-k v )/[l + (p f /p m )(l/X f -l)] 

kf +k m =1 

t( = 1 / 2 N f d f Jn ~Jk~f 


= 1/2 Jn/kf - 2_ 


d 


f 


3 5 =8^ 



k f =ti/4~0.785 


(a) Micromechanics and geometric relationships. 


Longitudinal modulus 


Transverse modulus 


Shear modulus 


Shear modulus 


Poisson’s ratio 


Poisson’s ratio 


Ei\ l kf E f 1 1 + k m E m 


' m 


V33 


1 - ^Jkf E f22 ) 


123 1 -ftf\'-G m /G f23 ) 

v m = kf Vfn +k m v m = v n 3 

v m = ' 


_ Ej22 i 


2 G 


£23 


3 



(b) Composite mechanical properties. 
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Longitudinal tension 

Longitudinal compression 
Fiber compression 

Delamination/ shear 
Microbuckling 

Transverse tension 
Transverse compression 
Intralaminar shear 
For voids 



Smc ~ 105^25 +2.5 S mT 


1 -kA 1 - 


SfllT ~ 1 — 


ifef - kf ) (l - E m j Ef22 ) S„ 


S (22 C ~ } ~ {jkf - kf ) (l - £■„,/ Ef 2 2 ) S n 


Sms ~ )(l-G m /Gfi 2 \ S„ 


S m *{l-[4k v /(l-k f )n] l,2 }s m 

(c) Composite uniaxial strengths, in-plane. 


2 



Void — 7 
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